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Abstract 

The prototypical Casimir effect arises when a scalar field is confined between parallel Dirich- 
let boundaries. We study corrections to this when the boundaries themselves have apertures 
and edges. We consider several geometries: a single plate with a slit in it, perpendicular 
plates separated by a gap, and two parallel plates, one of which has a long slit of large 
width, related to the case of one plate being semi-infinite. We develop a general formalism 
for studying such problems, based on the wavefunctional for the field in the gap between 
the plates. This formalism leads to a lower dimensional theory defined on the open regions 
of the plates or boundaries. The Casimir energy is then given in terms of the determinant of 
the nonlocal differential operator which defines the lower dimensional theory. We develop 
perturbative methods for computing these determinants. Our results are in good agree- 
ment with known results based on Monte Carlo simulations. The method is well suited to 
isolating the diffractive contributions to the Casimir energy. 



1 Introduction 



The Casimir effect, as originafly conceived, refers to the electromagnetic field in the presence 
of two infinite parallel conducting plates. The plates modify the boundary conditions on 
the field in a way which leads to a finite calculable shift in the ground state energy. Since 
the original work of Casimir [TJ, similar effects have been studied in a variety of different 
geometries, and a number of different calculational techniques have been developed. For a 
recent review see [2]. 

In the present work, we only consider scalar fields for simplicity. Our goal is to under- 
stand what happens when the boundaries themselves have edges. For instance, consider an 
infinite conducting plate with a hole in it. How does the size and shape of the hole modify 
the ground state energy? From the point of view of wave mechanics, the new feature is that 
the field can undergo diffraction as it passes through the hole. So our work could be viewed 
as the study of diffractive corrections to Casimir energies. 

An outline of this paper is as follows. We first develop a general formalism for studying 
edge effects, based on writing a lower-dimensional effective action for the field which lives in 
the hole. The total volume in which the field theory is defined is considered to be split into 
separate regions by boundaries, some of which have open regions to achieve the required 
geometry of surfaces. The lower-dimensional field theory is defined on the open regions of 
boundaries. This lower- dimensional action can be obtained by integrating out the scalar field 



in the bulk (section 2.1 ). We use this to study a single plate with a hole in it, and show that 



the Casimir energy can be expressed as the determinant of the nonlocal differential operator 



which defines the lower dimensional theory (section 2.2). Our effective action can also be 
described in terms of the wave functional of the field, projected onto the hole, providing a 
Hamiltonian type of interpretation where the normal to the surface takes on the role of time 



(section 2.3). In section [3] we specialize to a single plate with a long slit in it. We develop 
perturbative methods for computing the determinants, based on separating the operator 
into what could be considered as direct and diffractive contributions; mathematically these 
correspond to the "pole" (quasi- local) and "cut" (non-local) contributions in an integral 
representation. In section [4] these results are used to obtain the Casimir energy associated 
with two perpendicular plates separated by a gap. In section [5] we study two parallel plates 
separated by a gap, obtaining the diffractive contribution to the Casimir energy when one 
of the plates is semi-infinite. In both these cases, we compare the values obtained from our 
analytical calculation with the numerical calculations of the same available in the literature. 
The results from the two approaches are in good agreement. The appendices collect some 
mathematical results: the behavior of a field near a single plate with an edge (appendix 
|A| ) and heat kernels for the Laplacian with periodic and Dirichlet boundary conditions 
(appendix IbJ . 
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It is important to note that diffraction around solid obstacles does occur in many of 
the special geometries for which exact or close-to-exact results have been found. For the 
enormous amount of literature on such cases, we shall refer back to the reviews, except to 
point out that the multiple scattering techniques developed by a number of different groups 
[3] j [1] > [21 E] and the world line techniques of Gies et al [B] do incorporate diffraction around 
such solid objects. Nevertheless, the Casimir energy due to diffraction around edges of 
openings in boundaries have been calculated only in a few cases by world line techniques and 
Monte Carlo simulations [7]. The focus of our work is to develop an analytical understanding 
of such diffr active effects. 

2 An effective action for edge effects 

As a prototype for the sort of problem we will consider, take a free massless scalar field in d 
Euclidean dimensions. Imagine it propagates in two regions ("left" and "right") separated 
by a plate with a hole in it. Aside from the hole, we require that the field vanish everywhere 
on the boundary, while in the hole, we denote the fluctuating value of the field by <pQ. This 
is illustrated in Fig. [T] 

The basic idea is to write an effective action Sq for the fluctuations of <po . This effective 
action can be obtained in two different ways. From a path integral point of view it is a 
lower-dimensional effective action which arises from integrating out the scalar field in the 
bulk. From a Hamiltonian point of view, So is related to the wavefunctional of the field 
projected onto the hole. We develop the path integral approach first and return to the 
Hamiltonian approach in section [273| 



left 


right 








(|) = 



y = -b y = y = a 



Figure 1: The field propagates in two regions separated by a plate with a hole. At the 
location of the hole, indicated by the dotted line, we denote the fluctuating value of the 
field by 4>q. Elsewhere on the boundary, indicated by solid lines, we impose = 0. 
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2.1 Path integral approach 



We start with the Euclidean partition function 

Z = Jv^e'^'' 9 ^ 9 * (1) 
We fix the value of the field in the hole, </>|hoie = <t>o, and subsequently integrate over <fo- 

Z = Jvfo J V<f> e-S^kW* (2) 

0|hole = 4>0 

To perform the bulk path integral we set 

= <f>d + H (3) 

Here <p c \ is a solution to the classical equations of motion D^d = 0, subject to the boundary 
conditions 

0d = I 00 in h ° le (4) 
I elsewhere on boundary 

Since (f>o incorporates the boundary condition, 5(j) vanishes on all boundaries including the 
hole. The action for 5<j> can then be separated into left and right regions and the integration 
over this field can be done. This leads to 

Z = det- 1/2 (-a L )det- 1/2 (-D /? ) j Vfoe-i*** 8 **- 6 ** (5) 

where Or are Laplacians on the left and right regions. Given the boundary conditions 
on 54>, they act on functions that vanish everywhere on the boundary of the left and right 
regions (including the hole). 

To express <j) c \ in terms of (fio we introduce Gl and Gr, Green's functions on the left and 
right. They obey Dirichlet boundary conditions: they vanish everywhere on the boundary 
while in the bulk they obey ^lGl{x\x') = 5 d {x — x'), OrGr(x\x') = 5 d (x — x'). In terms 
of these Green's functions we have 

f j d d - l x' 4>o{x')n- d'G L (x\x') on left 

0d(x) = { (6) 
[ J d d l x' 4>o(x') n ■ d'Gn(x\x') on right 

Here n is an outward-pointing unit normal vector. Integrating by parts, the classical action 
in ([5]) is a surface term which can be evaluated with the help of ([6]). Putting this all together, 
we have 

Z = det- 1 / 2 (-D L )det- 1 / 2 (-a R ) yp0 o e- 5 ° (7) 
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where So = Sl + Sr with 



S, I d'' x x j d d - x x' ^ Q {x)M L {x\x') fotf) (8) 

S H - j d d - x x j d d - x x' M R {x\x') </> {x') (9) 



and 



M(x\x') = n- dn- d'G{x\x'), (10) 
appropriately for the left and right sides. 

The bulk determinants in ([7]) capture the Casimir energy that would be present if there 
was no hole. Corrections to this are given by a peculiar non-local field theory that lives 
on the hole separating the two regions; the fields </>o are nonzero only on the hole. We can 
write a mode expansion for the fields 4>o as ^o(^) = Ea c " u «( x ) where {u a (x)} constitute 
a complete set of modes for functions which are nonzero in the hole with the boundary 
condition that u a (x) — > as one approaches the edges of the hole. The action ^ takes the 
form 

Sl = ^ ^CaOLaP cp (11) 

where 

Laj3 = [ d A - x x d A ~ x x u a (x)M L {x\x) u p (x') (12) 



' hole 

with similar expressions for the right side of the partition. Because the mode functions 
u a (x) vanish outside the hole, this is essentially a projection of the operator M(x\x') to the 
hole. In other words, if we define an operator P which acts on functions / G L 2 (M. d ~ l ) by 

Pf(x) = l fix) ifxeh ° le (13) 
I otherwise 

then O = P M P. The functional integration now leads to 

Z = det " 1/2 (-Dl) det - l/2 {-U R ) det -^ 2 {O l + Or) (14) 

The explicit form of the operator M(x\x'), and its projected version O, will, in general, 
depend on the arrangement of plates and holes and boundaries. We can clarify the nature of 
this operator by constructing the Green's function with Dirichlet boundary conditions. For 
this, consider the right side of the box shown in Fig. [TJ We split the coordinates x = (x, y) 
into d — 1 coordinates x along the plate and a single transverse coordinate y. The plate is 
taken to be at y = 0. (Thus y is along the horizontal axis in the figure.) The right side 
of the box has length a along the y-direction, while we have lengths L±, L3, etc., along the 
other directions. The modes along the y-direction are 

MV) = \P- sinful (15) 



a \ a J 
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for n = 1, 2, etc. Similarly modes for the other directions take the form 



sm 



rriiiTXi 
Li 



with m.j = 1,2, etc. The Green's function for the left side can then be written as 



G(x|x') = -5> m (x)V4(x') 2 ~J2 



1 



. /niry 
a ^ p 2 + n 2 ir 2 /a 2 V a 



sm 



(16) 



(17) 



The summation over n can be carried out by complex integration (or other methods) to 
obtain 

G(x\x') = ~y~] - l -Vm(x)V'm( x/ ) coth(ap) [coshp(y - y') - coshp(y + y')] (18) 
lp 



where p = ^/p • p. This immediately leads to 



M R {x\x) = d y d y ,G{x\x) 

= ^^m(x)^(x') (pcoth(ap)) 



(19) 



Similarly, 



M L (x|x') = ^V m (x)^(x') (pcoth(6p)) 



(20) 



where 6 is the length (along y) of the left side of the box. For the parallel plate geometry, 
we are interested in the limit when b — > oo and L,- — > oo. The other cases we shall consider 



in this paper will also be special cases of the formulae (19) and (20). 



2.2 Single plate with a hole 

We can now go on to the projected version O of the operator M(x\x'). For this, we will 
first consider the example of a single plate with a hole in it. In this case, we are interested 
in a — > oo and b —> oo. Further, since Lj — > oo, we can approximate the sum over p by 
integration, to obtain 

M{x\x>) = J e^-^ VP^P (21) 

This result may also be obtained in a simpler way, without the full mode expansion, by 
noting that the standard Euclidean Green's function is 

G(x,y|x',y') ' 



(d - 2)vol(5 d - 1 ) (Jx - x'| 2 + (y - y r 



(d-2)/2 



6 



The Green's function appropriate to the plate geometry (Dirichlet boundary conditions at 
y = 0) can be constructed with the help of an image charge. 

G D (x,y\x',y') = G(x,y\x',y') - G(x, y\x', -y') 

The quantity we need is 

8 8 r I - 2 |x-x'| 2 -(ri-l> 2 

dydy,G D \ y=y , =0 - - vol( ^_ 1} (|x _ x/|2 + y2) (d + 2)/2 ^ 

where we have kept y — > + as a regulator. This is a quite singular-looking distribution: 
for x 7^ x' it approaches l/|x — x'| d , while at x = x' it diverges as —l/y d . 

To interpret this expression we return to the bulk equations of motion (V 2 + d y )<p = 0. 
A complete set of solutions is 

cl (x,y)=e ik - x e-^ for y > 

Notice that n ■ d <p c \ = —d y (j) c \ = k(p c i, so that we can identify 

n -d(j) c i\y = Q = \/-Vi (pd\ y =o 

Denoting the value of the field at y = by 4>q, and acting on ^ with n ■ d\ y= o, this implies 
that 

J d^x' cj, Q {x')d y d y ,G D {x,y\x' ,y')\ y=y , =G = T^V 2 o (x) (23) 



In other words, the distribution (22) is the square root of the Laplacian. Then the actions 
Q are 

S L = S R = f d d - l x ^qV^ (f>o (24) 
and the partition function Q is 



Z = det- 1 / 2 (-D L )det- 1 / 2 (-D i? ) J Vfa exp |- J d^xfayf^V* Q | 



(25) 



In this expression y 7 — V 2 refers to the square root of the Laplacian onl £i 1 , since that 
is what the arguments leading to ( 23 ) really establish^] But the path integral in ( 25 ) is over 



fields which vanish outside the hole. Denoting this qualification by the projection operator, 
the path integral can be evaluated to give 

Z = det - 1/2 (-D L ) det - 1 / 2 (-D i? ) det - 1 / 2 (P^V 2 P) 

The two bulk determinants can be absorbed by renormalizing the bulk cosmological constant 
and the plate tension]^] so the dependence on the size and shape of the hole is captured by 

Z = det ~ 1/2 (PV^V^P) (26) 
1 We are grateful to Alexios Polychronakos for discussions on this point. 

2 The stress tensor associated with an infinite Dirichlet plate is discussed in Birrell and Davies [5], section 
4.3. 
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2.3 Hamiltonian interpretation 



For an infinite plate with a hole our expression for the Casimir energy has a simple Hamil- 
tonian interpretation. Let us regard y as a Euclidean time coordinate. In an unbounded 
space the vacuum wavefunctional for the field is given by a Euclidean path integral over the 
region y > 0. 

<f>(x,y = 0) = 0o (x) 



Note that 0o is defined over the entire x plane. Following the logic in section 2.1 this means 
the vacuum wavefunctional is 

*oM = det- 1 / 2 (-n L )e- s ^ 



which, given (24), can be put in the familiar form [9| UOj 



*o[0o] = det- 1 /2(_ nL ) e -/^- 1 x^ov^v^ n 
So our expression Q for the partition function is really 

z = J P0 O m<M*o[M 

4>o = outside hole 
That is, to obtain the partition function we 



1. start with the vacuum state in the far past, at y = — oo 

2. evolve forward in time to y = 

3. impose a Dirichlet condition by only considering fields which vanish in the region 
outside the hole 

4. take the overlap with the vacuum state in the far future, evolved backwards in time 
to y = 

For a single plate with a hole our effective action is related to the vacuum wavefunctional 
by 

e -S»W=$*[^] vp o [0 o ]. 

A similar result holds in general, although with more complicated plate geometries one no 
longer has the standard vacuum wavefunctionals on the left and right. 

3 Plate with a slit 

For a single plate with a hole we have obtained a simple expression for the partition function, 

Z = det- 1 / 2 (P^V^P) 
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x = - a x = a 



Figure 2: In two Euclidean dimensions the slit geometry consists of two whiskers facing 
each other, separated by a distance 2a. 

where P is a projection operator onto the hole and \J — V 2 is the square root of the Laplacian 
on R^ -1 . To make further progress we now specialize to the case where the hole is a long 
slit of width 2a. The geometry is shown in Fig. [2j We first work in two dimensions, with a 
scalar field of mass fi. For a slit in higher dimensions, \x arises from Kaluza-Klein momentum 
along the transverse directions, so we can subsequently integrate over [i to obtain results 
appropriate to the dimension. 

In two dimensions the projection operator is 

p/0,0 = -?<*<• (27) 

I otherwise 

while V— V 2 can be defined by its spectral representation 

^V*e ikx = \k\e ikx k€R. 
However the operator we need to study is 



O = P^-V 2 + fi 2 P. 

This is not an easy operator to work with. In appendix [A] we diagonalize it for a slit of 
infinite widthj^] But for a slit of finite width we must resort to some sort of approximation 
scheme. 

To do this we note that by construction the field vanishes for x ^ a and x ^ —a. 
Moreover O respects a parity symmetry x —> —x. We therefore expect that we can expand 
the field in a complete set of odd- and even-parity functions which vanish for \x\ ^ a, namely 

| (— l) m -5= sin (ranxla) for — a ^ x ^ a , % 

^odd = ) y ' Va 'I ^ (2 8 ) 

otherwise 



„„„ f (— l) p+ 2 4= cos (p7rx/a) for — a ^ x ^ a , , 

otherwise 



The odd modes are labeled by n,m = 1,2,3,... while the even modes carry an index 
p,q = |,|, §,.... These modes are orthonormal; the factors of (— l) m and (— l) p+ 2 are 
inserted for later convenience. 



Meaning a single plate with an edge, described by Pf(x) = 



x < 
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These modes are eigenfunctions of the Laplacian with Dirichlet boundary conditions 
at x = a and x = —a. We will use them as a basis in which to diagonalize O = 
In this basis it turns out that O naturally splits into two pieces: a 
"pole" piece which is quasi-local and can be diagonalized, and a "cut" piece which is truly 
non-local. 



As a guide for the reader, in section 3.1 we consider the decomposition of O into its pole 



and cut contributions. In section 3.2 we set up the perturbation series and derive integral 



expressions for all higher terms in this expansion of the partition function. In section 3.3 
we integrate over the mass to find the ground state energy for a plate with a slit in four 
dimensions. 

3.1 Pole and cut contributions 



We begin by considering the odd-parity modes (28). They have a Fourier sine representation 



,odd _ 2 Va f°° Jlm rnr 



ib° = / dk , „ „ ^-^r sm(ka) sm(kx) 

tt ./q k 2 a 2 - n 2 7r 2 



Clearly = tp° , while for any function of the Laplacian we have 

F ( _ V 2 U° dd = — I" dk , 9 n ^ F(fc2 9 } 9 sm(ka) sin(kx) 
v ' tt Jo k 2 a 2 -n 2 TT 2 

It follows that the matrix elements in this basis are 
0£4 = < m |iV-V 2 + M 2 P\n) 



TT 



/ dk sin 2 (ka) m7T ^/k 2 +/j 1 2 



k 2 a 2 — n 2 TT 2 



Despite appearances, there are no singularities on the contour of integration: the would-be 
poles at k = ±rmr/a and irnr/a cancel against the zeroes of s'm(ka). 

Although O is not diagonal in this basis, the diagonal matrix elements are numerically 
much larger than the off-diagonal elements. There is a way of decomposing O which makes 
this manifest. First deform the integration contour slightly moving it just above the real k 
axis. Then write 

sin 2 (A;a) = - J {e 2ika + e~ 2ika - 2) 

For each term in this decomposition the integration contour can be deformed into the upper 
or lower half plane. One picks up a contribution if the integration contour crosses the poles 
(now real) at k = zizrmr/a or ±nn/a. One also gets a contribution when the contours 



4 One might question whether these modes provide a good basis in which to diagonalize O. This seems 
justified by the results of appendix [X] where we show that the exact eigenfunctions of O indeed go to zero 
as one approaches the edge of the slit (in fact they vanish as the square root of the distance from the edge) . 
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Figure 3: The contour for integrating over k. 

get wrapped around the cuts. The residues turn out to cancel unless m = n, so the pole 
contribution to the matrix element is diagonal. In fact 

odd,pole = I ^2^2+^2 (30) 

a 

The cut contribution to the matrix element is not diagonal. Rather we find 

«' cut = Hdy Vy^l (l - e-"^) 2 ™ 2 2 2 2 ™ 2 2 2 (31) 

it Ji m z TT z + fi z a z y z n z ir z + fi z a z y z 

Here y = lmk/fj, is an integration variable along the cut. 



Likewise the even-parity modes (29) have a Fourier cosine representation 



ipp Ven = [ dk 2 J 1 — cos(A;a) cos(/cx) 

7T Jq K CI — p 7T 

and the matrix elements of O in the even-parity sector are 

roo 



Deforming contours as before leads to the decomposition O pq = Opq le + 0™* where 

even,pole = I ^2+^2 ^ (32) 

' / <iy \/ y z — 1 (l + e ^ — „- 2 ■ p p „- — 2 „ ■ 2 „- p 

11 7r 7i p z 7r z + /j, z a z y z q z ir z + \x z a z y z 

Note the opposite sign in front of the exponential, due to the fact that the even-parity 
matrix elements involved cos 2 ka rather than sin 2 ka. 

Combining the even- and odd-parity matrix elements, note that pole can be identified 
with the operator yj — V 2 ^ + /i 2 , where V 2 ^ denotes the Laplacian with Dirichlet boundary 
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Figure 4: On the left, the geometry of interest. On the right, the geometry described by 
O pole , in which additional plates have been added to enforce Dirichlet boundary conditions 
at the edges of the hole. 

conditions at x = ±a. So although O pole is not a local differential operator, its square is 
local, and in this sense we will refer to O pole as being quasi-local. The cut contributions, 
on the other hand, make a truly nonlocal contribution to the operator O. 

From the physical point of view the decomposition into pole and cut contributions is 
natural because pole captures the geometrical optics effects of the hole, in which waves 
are directly transmitted from left to right, while O cut captures the diffractive effects. This 
follows from the observation made above, that O pole is related to an operator with Dirichlet 
boundary conditions at the edges of the hole. Such boundary conditions could be enforced 
by introducing additional plates as shown in Fig. |4| The additional plates prevent any 
diffraction from taking place, so diffractive effects are entirely encoded in O cut . 

3.2 Perturbation expansion 

Having decomposed the operator O = O pole + O cnt into direct and diffractive contributions, 
we wish to find a similar decomposition of the Casimir energy. This is straightforward. 
Expanding in powers of C cut , the partition function is 

-logZ = ^Tr log (0 polc + O cut ) 

= P ° 1C + ^Tr O p ^ c O cut - \Tr 0^ lc O cut O p ^O cut + ■■■ (33) 

The zeroth order term in this expansion gives the direct contribution to the energy, while 
the first and higher order terms give the diffractive contribution. 

Writing things in this way, the diffractive contribution to the energy is organized as a 
series expansion in powers of O cut . This expansion seems to be well behaved, even though 
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there is no small parameter in the problemj^] We give a speculative reason for this in the 
conclusions. But more prosaically the good behavior of the perturbation series will become 
evident from the explicit calculations we perform in the remainder of this paper, where we 
work up to 5 th order in O cut . For a graphical preview of the results see Fig. tal 



The lowest-order term in the perturbation series (33), the direct term, has a linear 
divergence and a subdominant logarithmic divergence, while all higher order terms, cor- 
responding to diffractive contributions, are logarithmically divergent. These logarithmic 
divergences are independent of a and can be eliminated by subtracting the a — > oo limit. 
This can be done either from the beginning, before the expansion in powers of O cut , or at 
the level of each term in the expansion. 

These subtractions can be interpreted as renormalizations of parameters corresponding 
to the plates and slits. Strictly speaking, in addition to the action for the fields, we have 
an action which describes the plates and slits in the given arrangement. This part of the 
action is generally of the form 

S = a A + a C + ••• (34) 

where A is the area of the plate, C is the length of the perimeters involved (for the plate and 
for any slits or holes in it). The coefficients a and a are the tensions for the plate and the 
edges of the slits. Being the coefficients of the area and perimeter terms, in the language of 
general relativity, they are the cosmological constants for the plate and for the boundaries. 
While these are calculable in terms of the material properties of the plates, at the level 
we are working, with the effects of the plates introduced as merely boundary conditions, 
they are free parameters. The partition function and the free energy we calculate are to 



be thought of as giving corrections to this action (34). The divergent terms we find can be 
absorbed as renormalizations of the parameters a, a. (In reality, at very short distances, 
the atomic structure of the plates become important and the divergent terms are rendered 
finite and calculable in terms of the interactions at that scale.) When there is a slit or hole 
in the plate, there is a part of the a A term missing and the renormalization of a appears 
in a way that depends on the dimensions of the hole or slit. This is because the term 
corresponding to the full area of the plates (ignoring holes and slits) is already subtracted 



out as explained at the end of subsection 2.2 and, so, the deficit is what is relevant for this 



part. The perimeter term should not depend on a, but only on the measure of the boundary 
(which is just two points for the one-dimensional slit) and can be identified easily by taking 
the large a limit. 



5 Note that the diffractive contribution to the energy does not in general have to be small compared to the 
direct contribution. Rather what we are claiming is that the diffractive contribution by itself has a useful 
series expansion in powers of O cut . 
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3.2.1 The direct contribution (lowest order) 



We now proceed to study the various terms in (33). Combining the even- and odd-parity 
contributions, at leading order we have 

-logZW = \ logV(W«) 2 + /" 2 

= l -TY\og(-V 2 D + n 2 ) (35) 

where again V|j is the Laplacian in the slit with Dirichlet boundary conditions. This can 
be computed using the heat kernel methods described in appendix [B) There is a linear di- 
vergence, proportional to a, which renormalizes the cosmological constant in the slit. There 
is also a log divergence, independent of a, which renormalizes the boundary cosmological 
constant (i.e. the tension associated with the edges of the slit). After these divergences 
are removed one is left with a finite result which vanishes exponentially as a — > oo. In the 
notation of appendix [Bj 



00 ds 



2 

sjj, 2 \ A — 4a 2 n 2 /s 



n=l 

ilog(l- e - 4 ^) (36) 



3.2.2 The diffractive contributions 

For the first order term it is useful to separate the contributions from the odd and even 
parity terms as 

- log ZW = - log Z« - log Z« n (37) 

where 

° dd VT A ^ ^/ r 2^2 + M 2 a 2 ( r 2^2 + M 2 a 2 y 2)2 

-logZ« n = -^-/ dyV^T J] \ ^ — — (39) 

For ease of presentation of these and higher order results, we define 
T(pa,y,z) = V 

7T *— * 



7T y 7 r 2 7r 2 + ^ 2 a 2 (r 2 7r 2 + fi 2 a 2 y 2 ) (r 2 ir 2 + fi 2 a 2 z 2 ) 

roc r 2 7r 2 

^ J-oo dX ( r2 K 2 + /U 2 a 2 + A 2 ) (r 2 7r 2 + p?a 2 y 2 ) (r 2 n 2 + fi 2 a 2 z 2 ) ^ 
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By resolving the integrand into partial fractions, the summation can be done using 

1 



Y — 



vr 2 + A 2 2A 2 



[AcothA - 1] 



(41) 



We can then write T as 



CO 



T(fia,y,z) = -a / d\ 
^ Jo 



+- 



VX 2 + 1 coth(/WA 2 + 1) 
~(A 2 + l-y 2 )(A 2 + l-z 2 ) 

y coth(//ay) z coth(^az) 



+ 



(y 2 -z 2 )(A 2 + l-y 2 ) (z 2 -y 2 )(A 2 + l-z 2 ) 
The limit of z — > y is seen to be 

1 d rco 



(42) 



T(fJ,a,y,y) 



2jr 2 y dy 



dX 



-y/X 2 + 1 coth(/xaV A 2 + 1) + y coth(fiay) 
A 2 + 1 - y 2 



(43) 



The a — > oo limit of these expressions is logarithmically divergent. The renormalized 
contribution is obtained after subtraction of this divergence as 



- log Z {l ) A 

o odd,ren 



/oo 
dyVv 2 ~ 1 [(1 " e~ 2/ ^) r(^o, y, y) - T(jia -> oo, y, y)] (44) 



Similarly, for the even-parity contribution, we define 



2„2 



S(jm,y,z) 



pL a 



E 



r 2 7r 2 



_! 3 ^/ r 2 7r 2 + fi 2 a 2 (r 2 7r 2 + fi 2 a 2 y 2 ) (r 2 n 2 + fj, 2 a 2 z 2 ) 



(2yua) 



7T- 



E 

Z=odd 



dA 



(Z% 2 + (2^a) 2 + A 2 )(Pvr 2 + (2//a) 2 y 2 )(/ 2 vr 2 + (2//a) 2 z 2 ) 

(45) 



Again, by use of partial fractions and (|41|), we can write this as 

1 /"°0 



S(fia,y,z) = -g / dA 
^ JO 



VA 2 + 1 [coth(2/WA 2 + 1) - csch(2^a v / A 2 + 1) 
(A 2 + l-y 2 )(A 2 + l-z 2 ) 



y[coth(2/iay) — csch(2/xay)] z[coth(2^az) — csch(2fiaz)} 



(y 2 - z 2 )(A 2 + 1 - y 2 ) 



{z 2 - y 2 )(A 2 + 1 - z 2 ) 



(46) 



S(na,y,y) 



1 d 
2?r 2 y dy y 



+ 



VA 2 + 1 [coth(2/WA 2 + 1) - csch(2/WA 2 + 1)] 
A 2 + 1 - y 2 

y[coth(2^ay) — csch(2/xay)] 
A 2 + 1 - y 2 



(47) 
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Figure 5: Free energy vs. fia for a two dimensional slit. Lower curve is upper 
dashed curve is - log and upper dotted curve is -logZ( 2 ). 



The renormalized expression for the even-parity contribution is then 

/oo 
dy Vv 2 ~ 1 [(1 + e- 2 » ay ) S(»a, y, y) - S(fxa -> oo, y, y)] (48) 

The higher order terms can also be written down easily in terms of T(fia, y, z) and 
S(fia, y, z) as 



1o S Z odd 



in— 1 poo 



J) 



JJ(1 _ e - 2 ^)T(/xa) * T(/ia) * • • • * T(/ia) 



— (/ixa — >• oo) 



(49) 



log Z (ri) 



in— 1 /-oo 



n 



/oo "■ , 



+ e- 2 ^)5(/ia) * S(/xa) * • • • * 5(/xa) 



-(/ia — > oo) 



where 



T * T * ■ ■ ■ * T = T(pa, yi,y 2 )T(fia, y 2 , 2/3) •• • T(fia, y m yi), 
with a similar expression for the S"s. 



(50) 



(51) 



The integrals involved in these formulae can be computed numerically as a function of 
[ia. The direct term and the first two diffractive contributions are shown in Fig. [5] Notice 
that the second order diffractive term is much smaller than the first order term, consistent 



with our expectation of the usefulness of the expansion ( 33 ) . 
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3.3 Slits in 4 dimensions 



We can now extend our results to the physical setting of four dimensions by introducing two 
more dimensions: a periodic Euclidean time dimension of size f3 (representing the inverse 
temperature), and a space dimension of size L (representing the length measured along the 
edge of the slit)j^] For large j3 and L the 4-dimensional partition function is an integral, 

- log Z M = (3Lj (- log Z 2d ) (52) 

where we are interpreting the momentum in the extra dimensions as providing a Kaluza- 
Klein mass. This means the energy per unit length for a slit in four dimensions is 

E r ^(-lo g Z 2d ). (53) 



L /n 2tt 



The direct contribution is thus given, using (36), by 

C(3)L L 



E< °> = - 5 ( - 2 - 99 * 10 ~ 3) (54) 



The diffractive contributions, obtained by integrating — logZ( n ' from (44), (48), and (49), 
(50) are 

= L (2.15 x 1(T 3 ), £ (2) = 4 (0.14 x 1(T 3 ) (55) 
a z a z 

The total value of the energy, to this order, is —0.7 x 10~ 3 (L/a 2 ). The 1/a 2 dependence of 

these results is, of course, fixed by dimensional analysis. 

We may also note that the energy for a slit in arbitrary number of dimensions can be 
obtained by extending the integration over //in (52) to higher dimensions. 



4 Perpendicular plates 

In this section we study the ground state energy for two perpendicular plates separated by 
a distance a. The geometry is shown in the left panel of Fig. [6j 

Again our basic approach is to write a lower-dimensional effective action in the gap 
between the plates, indicated by the dotted line in the figure. Fortunately this turns out to 
require very little effort. The odd-parity modes in a slit that we studied in section [3] vanish 
at x = 0. Thus they are appropriate for the case of there being a plate perpendicular to the 
slit as in the right panel of Fig. [6} We are interested basically in one side of this geometry. 

6 While we use periodic boundary condition for the time direction, we will retain Dirichlet conditions for 
the spatial directions. In the limit of large L, we can replace summations over momenta along this spatial 
direction by integration. The distinction between Dirichlet conditions and periodic boundary conditions will 
not matter as L — > oo. 
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x = a 


- a 


a 



Figure 6: On the left, perpendicular plates separated by a distance a. On the right, the 
geometry described by the odd modes in a slit. 



Notice however that the modes relevant for this case, namely, for $S x $J a are in one- 
to-one correspondence with the modes relevant for the slit —a^x^a. The eigenvalues 
are also the same. The determinant is then given by the result for the odd-parity modes of 
the slit. So to obtain the Casimir of the two perpendicular plates, we can take the result 
for — log Z2d discussed in the previous section, but restricted to the odd-parity modes, and 
then integrate over fi for the appropriate number of transverse dimensions Q Carrying out 
the integrations in four dimensions, we find 



Ef = 


L 




4 1 ) = 


L 


a" 2 


E f = 


L 


~a? 


E f = 


L 





(-11.96 x 1(T 3 ) 
( 5.01 x 1(T 3 ), 
( 0.16 x 10~ 3 ), 
( 0.01 x 10~ 3 ) 



E 



(2) 



L 

ft- 



0.66 x 10~ 



E 



(i) 



L 

n - 



0.05 x 10" 3 ) 



(56) 



The total value for the E± up to this order is —6.07(2) x 10 _3 (L/a 2 ). The terms in (56) 
have been evaluated using Mathematica. The n-th order term involves 2n + 1 integrals. 
As the number of integrals increases, the precision of the answers is lowered. We used 
several integration methods suitable for multi-dimensional integrals. Comparing results 
from different integration methods we estimate the error in our final answer for the total 
value to be within 0.02 x 10~ 3 . 

The Casimir energy for two perpendicular plates separated by a gap has been numerically 
investigated by Gies and Klingmuller [7]. Their calculation is done by considering a path 
integral representation for the propagator. When the two plates are present, all paths which 
touch both plates must be considered as an overcounting of paths and must be removed 
from the sum over paths. This process, in a Monte Carlo evaluation of the path integral, 
then leads to corrections to the pure vacuum result and gives the Casimir energy. Their 



7 The result for the full geometry of the right panel will require independent modes for the left and right 
sides of the vertical plate. We must use modes sin(n7r:r/a), sin(m7rs:/a) with m,n being independently 
chosen integers. This will lead to a doubling of our results for that geometry. 
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final result for two perpendicular plates is given as 

L 



Ei 



(-6.00(2) x 10~ d ) 



(57) 



Clearly our result is in very good agreement with the above value calculated in [7J. 



5 Parallel plates 

Consider an infinite plate with a hole in it, parallel to a second infinite plate with no hole. 
Let the separation distance between the plates be a. The geometry is shown in Fig. [JJ 

It is straightforward to study this situation along the lines of section [2| We are interested 
in keeping a finite but taking b — )■ oo. In this case 



M L (x\x') + M R (x\x') = V>m(x)V4(x') P (1 + coth(ap)) 



(58) 



Since V'm(x) form a complete set of states, this is basically the operator \J — V 2 (l+cothav / — V 2 ). 
Once again, an alternative way to arrive at the above equation is the following. A complete 
set of solutions to the bulk equations of motion in the region between the plates is 

<Aci(x,y) = e ik - x (Ae- k y + Be k y 

where the Dirichlet boundary conditions at y = a fix 

B = -Ae~ 2ka 



Note that 



n ■ 



M\y=o 



y<Pcl\y=0 



k(A-B)e 



ikx 



k 



A-B 
A + B 



M\y=0 



x = W 
x = 



y = 



y = a 



Figure 7: Parallel plate geometry. The z axis is into the page. 
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This means that we can identify 



n -^ dl - = tanh(a^r lls/ =° 
The actions ([8]), ^ are then given by 

S L = j rf^-'x^oV^V 1 ^ (59) 



This can also be interpreted in the Hamiltonian language of section 2.3 The wavefunctional 
on the left has the standard vacuum form 

* L [M = det-^ 2 (-n L ) exp (- j rf^xi^oV^foJ 

while the presence of the second plate modifies the wavefunctional on the right to 

V R [<fo\ = det - 1/2 {-U R ) exp (- J d*- 1 * ^oV /= V2 coth(aV-V*)<h) 

From either perspective, the partition function ([7]) is 
Z = det- 1 / 2 (-D L )det- 1 / 2 (-D iJ )det- 1 / 2 [p(7^V2 + V^V2coth(ay^V2))pl , (61) 



in agreement with (14) and (58). Just to be clear: the first determinant is computed in 
the region y < with a Dirichlet boundary condition at y = 0. The second determinant 
is computed in the region < y < a with Dirichlet boundary conditions at y = and 
y = a. In the third determinant P is a projection operator onto the hole and y/— V 2 is the 
Laplacian on 

The first determinant, det -1 / 2 (— □/,), can be renormalized away. The second determi- 
nant, det- 1 / 2 (-D jR ), gives the Casimir energy that two parallel plates would have if there 
was no hole. In 3+1 dimensions this Casimir energy, per unit area, is given by the standard 
result 

E 7T 2 

A ~ "1440^3 ( 62 ) 
But in our case, for the sake of comparison with numerical results, it is important to specify 
boundary conditions at infinity (meaning on the walls of the box shown in Fig. 0. In 
that figure, if we increase the size of the hole until it reaches the walls of the box, the 
condition that vanishes at the edge of the hole is carried over to a Dirichlet condition 
on the walls of the box. So increasing the size of the hole is consistent provided we use 
Dirichlet boundary conditions at infinity. That is, det _1 / 2 (— should be computed in 
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a box of size L\ x a x L3, where L%, L3 — > 00 are the lengths in the x, z directions. This 
leads to subdominant terms in the Casimir energy, namely 



E 



-J^_L X L Z + -®-(2Li + 2L 3 ) + 
1440a 3 d 64vra 2V d; 



(63) 



Other choices of boundary conditions at infinity are possible. For instance, if we use periodic 
boundary conditions in the z direction, the term £(3)Li/327ra 2 is absent and we'd have 



E 



7T 



1440a 3 32vra 2 
For simplicity this is the case we will treat in the following. 



(64) 



Going back to (61 ), we note that all the dependence on the size and shape of the hole is 
captured by the third determinant, which we now proceed to study. To keep the discussion 
simple we take the hole to be a long slit of width W and length L = L3. We expand the 
field in the slit in modes analogous to ( 28 ) , ( 29 ) , 

~2 



The slit is located at < x < W so that n 



sin{niTx/W) 



1,2,3,.... 



(65) 



As in section 3.1, the operator we are interested in, namely, 

O = P V^V 2 (l + coth(aV-V 2 )) P (66) 

can be decomposed into pole and cut contributions in this basis. The free energy F = 
Tr logO can then be expanded in powers of O cut . The matrix elements of the projected 
operator in (66) are given by 

4 rW 
W 



dxdx 1 







dP e i P {x-x') f $n( mrx /W)Bhi(mira//W) f(p) 



VP 2 + A* 2 



1 



1 — exp(—2a\/p 2 + n 2 ) 



(67) 



To evaluate the cut-terms arising from the square root factors, it is useful to write an 
integral representation for f(p), namely, 



M = - 

7T 



da P 2 + ^ 2 

-00 (a 2 +p 2 + M 2 ) 



1 



1 — exp(2am — e) 



(68) 



The exponent e pushes the poles at a = mr/a to the lower half-plane. We can evaluate the a- 



integral by completing the contour in the upper half-plane; only the pole at a = i\/p 2 + /J? 
contributes and the equivalence with (67) can be easily verified. Carrying out the integration 
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over x, x', we then find 

n 4 f dp 

^mn 

H(q,q') 



W J 2vr 
1 



f(p)H(q,q') 



I _|_ e i(q-q')W _ e i(p-q)W _ e -i(p-q')W 



(p - q + ie)(p - q' - ie) 

-(g -> q' -> g') - (g g, g' -> -g') + (g 



-g,g 



(69) 



where g = nir/W, (/' = rmr/W. The integration over p can now be done. For a term with 
1 + e*(9~9')W — e *(P — we nee £ Q c i ose ^he contour in the upper half-plane, for a term 
with e~ i ^ p ~ q ' W , we need to close in the lower half-plane. There will be pole contributions 
from the denominators p — q + ie and p — q' — ie. These are identical to what we named 
the pole terms in O mn . There will also be terms from the poles of f(p). The latter will 
correspond to the cut-terms we are seeking. The evaluation of the p-integral then leads to 

n _ /ripole 



CO wit 



^)pole 



2u{q) 



cut 



A{a,q,q) 
where uj(q) 



irW 

oo 



'mn = 2/(g) d m 
[l + (-l) m+n }A(a,q,q') 



\ _ e —2au)(q) 

Aqq' 



(70) 



dX 



/(A) 



+ 



/(g) 



+ 



Z(g') 



(g 2 -A 2 )(g /2 -A 2 ) {q> 2 - q 2 ){\ 2 - q 2 ) {q 



,/2 



)(A 2 



y/q 2 + H 2 . 



5.1 The direct contribution (lowest order) 

The direct contribution to the free energy is given by 

f3F = ^TrlogO pole (71) 
= I T rlog(2^V 2 ;)-^Trlog(l-e- 2a v /Z Vi) (72) 

where V 2 ^ is the Laplacian in the slit with Dirichlet boundary conditions at x = and 
x = W . The first term has UV divergences but is independent of a. In fact it is just the 



lowest order energy in a single slit which we studied in section 3.3 There we found that 
the energy per unit length for a slit of width W 

(73 ) 



^Previous versions of this paper had an overall sign error in 

9 One can obtain this directly, as E/L = (1//3L) | Tr log(— V|>), using the results in appendix B 
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The dependence on the separation between plates is captured by the second term in (72) 
which is finite in the UV. Including both terms, we have the finite (renormalized) energy 
per unit length 



^ = C(3) _ 1 f A / _ 2n ^HW ) (7D 

L 32irW 2 2 J (2vr) 2 2^ g V J 1 ' 



This expression can be studied in various limits. As W — > the second term makes an 
exponentially small correction, and we have 

— ~ . 1 r -2na/W (7 r\ 

L~ 32ttW 2 8aW 1 1 

On the other hand as W — > oo the second term dominates. To study it in this limit we use 
the Euler-Maclaurin summation formula, 

22f{w) =w J ^/(^)- 2 /(o)-E (2m)! w2m ^ co) 

This leads to 



n=l m= 



^ (2m)! W 2m ~ l 

E tt 2 W C(3) 



L 1440a 3 32vra 2 



+ (77) 



Adding this result to the bulk contribution (64), with L\ and L% = L taken to be large, we 
find, for the lowest order or direct contribution to the energy, 

^ 2 (L X - W) 



L 1440a 3 



+ ••• (78) 



The terms proportional to 1/a 2 cancel out. 10 Also the usual Casimir energy per unit area 



(62) in the region corresponding to the slit is canceled out and only the facing area of the 



two plates L(Li — W) appears in E^°\ 

5.2 First diffractive contribution (first order) 

The diffractive contribution to the 2d free energy arises from the expansion 

- [log Z - log Z(°>] = \ Tr log fs mn + ^ [1 + (-1)™+™] A(a, q, q')) (79) 

We can easily work out the higher order terms from this. In the case when W is large, 

2L f°° 1 — e~ 2cUJ (.i) 

-logZW = — / dqq 2 — A(a,q,q) 

K l Jo 



id 



This cancellation depends on the boundary conditions at infinity used in ( 64 1 , namely Dirichlet in x and 



periodic in z. With a different choice of boundary conditions at infinity the 1/a terms would not cancel 
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The n-th order term is given by 



™ j poo 

J _■ 



2 n L r°°-fr, 2 A («) A («) A (°) 



/ / / 

A(a) A (a) A(a,gi,g 2 ) A(a, g 2 , 93) A(a,g n ,gi) 



* • • • * 



/ / f(qi) f(Q2) f(q n 



(81) 



As in the case of the slit and the perpendicular plates, the renormalized expressions are 
obtained by subtracting the a — > oo limit. 



on j poo n 

iog^ = (-ir +1 ^/ o U d ^ 



A(o) A(a) A(a) , , 

— ^ * — ^ * • • • * — Y~ - (a^oc) 



(82) 

The energy for the case of four dimensions can now be obtained by integration over /i, 

E = J ^(-logZ rcn ) (83) 
Evaluating the integrals numerically, we find, for the first few orders, 

= ^ ( 5-54 x 1(T 3 ), = ^ ( 0.80 x lO" 3 ), 

£ (3) = 4 ( 0.19 x 10~ 3 ), £ (4) = ( 0.05 x 10~ 3 ), 
or or 

E^ = — ( 0.01 x 10~ 3 ) (84) 
a 1 

The value for the diffractive contribution to the energy, up to this order, is E dlSr = 6.59 x 
10~ 3 (L/a 2 ). To obtain the total energy associated with the edges of the slit we combine the 



direct contribution which appears in (77) with the diffractive contribution computed here, 
to fined] 

E siit = _ C(3)L + EdiSr ^ _ 5 3? x 10 _ 3 L 

327TGT or 

This result is for a slit of finite, although large, width. There is no direct comparison to 
other methods of calculation available. However, the case of two parallel plates, one of which 
is semi-infinite, provides a point of comparison. The Casimir energy for this geometry has 
been numerically investigated by Gies and Klingmuller [7], by the world-line method of 
subtracting out the paths which touch both plates. Their final result for two parallel plates, 
one of which is semi-infinite, is given as 



I or 

where 7 = 0.00523(2). 



11 Sin ce w e are interested in the energy associated with the edges of the slit, we leave out the usual Casimir 
energy (781 associated with the facing area of the two plates. We also leave out the £(3)Z/3/32-7ra 2 term in 



(64 1 since, rather than being associated with the slit, it is associated with the corners which appear at 



infinity when a Dirichlet condition at infinity in the x direction is imposed. 
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We have a slit of finite, although large, width. Thus there are two edges to the slit, 
each of length L which must be considered. If one were to remove paths from a sum-over- 
paths formula for the propagator, all paths which touch both edges must be removed. The 
calculation in [7J is for a semi-infinite plate to begin with and hence paths which touch on 
the edge which is far away from x = are not removed. Thus our result must be divided 
by 2 for the edge terms to get a proper comparison, and we find 

E eAze = _lL w . th 5 37x 1 q-3 (87) 

2 a z 

This value for 7 is in good agreement with the result of Gies and Klingmiiller [7J . The inclu- 
sion of higher order diffractive terms would further decrease our value for 7 and presumably 
bring it closer to the world-line result of |7j. 

An exact calculation of the Casimir energy of a parabolic cylinder next to an infinite 
plate has been done by Graham et al [11] . A particular limit of this gives the result for the 
case we are studying, namely, a semi- infinite plate next to a parallel plate. The result in 
|llj is I7 = 0.0025, again in agreement with our result and with [7j. 

6 Summary 

We have developed a method for calculating Casimir energies, including diffractive con- 
tributions which can arise from apertures on plates and other boundary elements of the 
geometry. This involves the functional integration over a lower dimensional field theory de- 
fined on such apertures. The relevant kinetic operator has an interesting structure. In the 
simplest cases it is of the form V— V 2 , similar to what occurs in the wave functional of the 
field, but there are important modifications based on the geometry of the situation. In all 
cases, the operator acts on functions which have support only on the apertures. The matrix 
elements of the operator allow a clean separation of diffractive contributions from direct 
(or ray optics) contributions. To evaluate the relevant functional integrals we expanded in 
powers of the diffractive contribution. This seems to be a good approximation even though 
there is no explicit small parameter in the problem. 

In this paper we focused on the Casimir energy for some special cases: a single slit, two 
parallel plates, one of which has a long slit in it, and two perpendicular plates separated 
by a gap. In the latter two cases numerical calculations based on world line methods have 
been performed. Our results can be compared, and in both cases the agreement is quite 
good. But the method we have developed is quite general and can be applied to a variety 
of different geometries. For instance it can be easily generalized to arbitrary dimensions. 
In fact working in d dimensions might justify the perturbation series, as an expansion in 
powers of 1/d. The method could also be extended to include finite temperature effects, 
which would allow a comparison with the results of [14J. 
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A Exact mode functions on an interval and half- line 



In the formalism we have developed, effects associated with a hole are captured by non-local 
differential operators of the form PFP where P is a projection operator onto the hole and 
F is some function of the Laplacian on R ' -1 . Diagonalizing such operators is in general 
difficult and we were forced to resort to perturbation theory. However one can find the 
exact eigenfunctions numerically. Also in some cases an analytical treatment is possible, 
using results obtained long ago by Malyuzhinets for scattering from a wedge |12} IT3] . Here 
we collect some of these results. Besides illustrating the non-perturbative features of the 
problem, our motivation is to provide evidence that a perturbative treatment should be 
reliable. 



First consider the operator PFP where F = \J— d 2 /dx 2 and P is a projection operator 
on the unit interval [0, 1]. One can study this numerically, starting from a finite difference 
approximation to the Laplacian in position space. 



(F 2 



{F 2 )i,i+i 



(F 2 



In this basis 





( 





\ 


p = 





1 












0/ 



One can take the square root of (88) numerically and diagonalize PFP. A typical eigen- 
function is shown in Fig. [8} The exact eigenfunction clearly vanishes at the edges of the 
interval and approaches a plane wave in the middle. It is well- approximated by the Dirich- 
let modes sin(n7r:z;) that we used as the basis for our perturbation series. Indeed the only 
significant difference between the exact and perturbative modes is that the exact modes go 
to zero more steeply at the edges of the interval. 



To study this edge behavior we consider F = 
operator onto a half-line. 

/ fix) 
) 



cP 
dx 2 



+ fj? and take P to be a projection 



Pf(x) 



x > 
otherwise 



(89) 
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This describes a semi-infinite "plate" embedded in two Euclidean spacetime dimensions. As 
this geometry has no adjustable parameters, up to renormalization the energy of such a plate 
vanishes, and in this sense there is nothing interesting to calculate. Instead our motivation 
for studying this geometry is that we will be able to diagonalize PFP analytically and 
show that the exact eigenfunctions have y/x behavior as x — > 0. We will show this in two 
different ways: first using Laplace transforms, then by solving the wave equation following 
Malyuzhinets. 

A.l Laplace transform 



The eigenvalue problem we wish to solve is 

PFP(j)(x) = X4>{x) 



One way to solve (90) is to find a function 4>{x) such that 



(90) 



1. 4>(x) has support for x > 0, so that Pcf) = eft 

2. (F - \)4> has support for x < 0, so that P(F - X)(j) = 

Suppose we represent 4>(x) using an inverse Laplace transform. 

p<0 ° ds 



4>(x) 



2-rti 



e sx <t>(s) 



(91) 



If (j> is analytic for Res > 0, then <p(x) will vanish for x < 0. Likewise if <f> is analytic for 
Res < 0, then <j){x) will vanish for x > 0. So we need to find a function 0(s) such that 




Figure 8: The 5 mode in a slit. The red solid curve is the exact eigenfunction. The black 
dashed curve is sin(57rx). 
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1. <j>(s) is analytic for Res > 

2. (y/ (j? — s 2 — \)4>(s) is analytic for Res < 



It's convenient to work on the covering space of the cut s plane, making a change of variables 



s = i^isinhz. We also set A = \i cosh £ 12 Then (91) becomes 

Mx) = — r dz coshze^ sinh2 0(z) . 

27r J-oo 

If 4>(z) is analytic in the strip — n < Imz < and satisfies 

(j)(z) = (p(—z — in) 

then we can replace 

dz ^n( + +/ +/ ) dz 

3 " \J — oo J oo Joo—in J — oo—iir J 

and (j){x) will vanish for x < 0. (The conditions on (j)(z) correspond to the requirement that 
(f)(s) is analytic and single-valued for Res > 0.) Likewise if (f)(z) is analytic in the strip 
< Imz < n and satisfies 

4>(z) = (p(—z + in) 

(corresponding to the requirement that <j>(s) is analytic and single- valued for Re s < 0) then 
4>(x) will vanish for x > 0. So corresponding to the conditions on 4>(s), we need a function 
4>(z) such that 

1. 4>(z) is analytic for — n < Imz < and satisfies 

fa) = 4>(-z - in) (92) 

2. (coshz — cosh^)0(z) is analytic for < Imz < n and satisfies 

(coshz — cosh£)0(z) = — (coshz + cosh£)0(— z + in) (93) 

A solution to this system of equations was obtained by Malyuzhinets |12} [13] . Define 

*7r(^) = Vtt( - iz + y +^)V , 7 r( - + — - - iz - - + i()^7r{ - iz - - - i£) 

where 

- exp--^ S cosh(t7r/2) sinh(2at) ' 



12 



The 5-function normalizeable spectrum of PFP is A £ [/i, oo) corresponding to < £ < oo. 
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a 

D 



Figure 9: The wedge geometry. We impose Dirichlet boundary conditions on the plate at 
9 = 0. To diagonalize the operator PFP we impose Robin boundary conditions on the hole 
at 9 = a. 

Then the solution is 



(sinh z — ie) 2 — sinh £ 



To see this, note that ^f n (z) was constructed to satisfy (92), (93) all by itself. Moreover the 
prefactor l/(sinh 2 2: — sinh 2 £) is invariant under z — > —z ± in, so 4>(z) also satisfies (92), 



(93). It only remains to check the analyticity conditions. In the strip —n < Imz < n one 
can show that ^ n (z) has no poles. It does, however, have zeroes at z = ±£ + in. In the 
same strip the prefactor has poles at z = ±£ + ie and z = ±£ + in — ie, but the latter poles 
cancel against the zeroes of ty n . So in the strip <p has poles at z = ±£ + ie. That is, 4>(z) 
is analytic for — n < Imz < and, when multiplied by coshz — cosh£, it becomes analytic 
for < Imz < n. 



One comment on this solution is in order. Up to a normalization the prefactor in (94) is 
the Laplace transform of 9{x) sin kx, where k = M sinh£ - exactly the modes with Dirichlet 
boundary conditions that were the starting point for our perturbation theory. So in the 
case of a semi-infinite plate, diffractive corrections to the perturbative modes are given by 
the Malyuzhinets function ^>- K {z). 

A. 2 Wave equation in a wedge 

Another approach to diagonalizing PFP more closely makes contact with the original work 
of Malyuzhinets. Consider a semi-infinite plate in two Euclidean dimensions, and let's 
study the wavefunctional for the field on a 'hole' which is a half-line making an angle a 
with respect to the plate. In general the surface actions appearing in Q, ^ can be written 

as 

S = / -(j> c \n ■ d $ c \ 

J hole 1 

where (p c \ is a solution to (— □ + fj?)<p c \ = with the boundary conditions <p c \ = at 9 = 0, 
<t>d = (fro &t 9 = a. Rather than specify the value of <fio, suppose we impose n ■ d(j) c \ = Xcj) c \ 
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at 6 = a. That is, suppose we solve the system 

(-□ + /i 2 )0 cl = 

c i = at 6 = (Dirichlet boundary condition) (95) 
n ■ d(j> c \ = \<t> c \ at 9 = a (Robin boundary condition) 

To make contact with the problem of diagonalizing PFP, note that when a = n we have 
another expression for the surface action, namely 

hole 1 

So when a = tt, the value of the field along the Robin boundary <j) c i(r,9 = a) is a solution 
to PFP(f> = \cj>. 



Fortunately ( 95 ) is exactly the problem studied by Malyuzhinets [12|, [T3] . For general 



a, and denoting A = ^cosh£, the solution is 

Mr,6)= [ ^e- mrcos ^f(z)g(z) 

where 

f(z) = il) a (z + a + i£)i>a(z + a - i£,)ip a (z - a + i£)ip a (z - a - it) 
is the function introduced by Malyuzhinets and 



sin 2 (z^z) — sin 2 (z//3) 

is chosen to obtain the correct asymptotic behavior at large r. Here v = n/2a and /3 = 
| — a + i£. The contour 7+ starts at z = 2ir + ioo, descends towards the real axis, moves to 
2tt to the left while staying above any singularities of the integrand, and returns to +ioo, 
while 7_ is the mirror image of 7+ under z — > — z. 

One can extract the asymptotic behavior of the solution from the contour integral rep- 
resentation |12[ 113] . Along the Robin boundary <f> c \ has plane- wave behavior at large r, 
while near the origin it has power- law behavior c i( r ) ~ r n / 2a . Setting a = tt, this means 
eigenfunctions of PFP have plane-wave behavior far from the edge, while near the edge 
they vanish like \fr as r — > 0. 

B Heat kernels 

For a free scalar field of mass m in Euclidean space, with some number of periodic dimensions 
and some number of Dirichlet directions, the partition function is 

Z = e- /3F = det- 1 / 2 (-D + m 2 ). 
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Here /3 is just the periodicity around some "Euclidean time" direction. We can represent 



1 



Tr log(-D + m 2 ) 



ds 



Tre~ s(_D+m2) 



1 

~2j e 2 s 

where e — > serves as a UV regulator. To compute the (trace of the) heat kernel 



(96) 



K(s) = Tre' 



we use the fact that □ = gj^ where the eigenvalues of df are 



Xi periodic with period Li 
Xj Dirichlet with size L; 



4> eigenvalues 

eigenvalues n£N 
This means that the heat kernel factorizes, K(s) = Y\i^i^ where 



Xi periodic 
Xi Dirichlet 



Note that 



K i = Kp(s,L l ) 
Ki = K D (s,Li 

1 



E e 

nSZ 



-s(2nn/Li) 2 



-s(nn/Li) 2 



K D (s,L) = -(K P (s,2L)-l) 



By Poisson resummation 

K P {s,L) 



L 



'Airs 



-L 2 n 2 /4s 



0, 



L 2 



(97) 
(98) 



The expressions (97), (|98j) make it clear that as s — > we have 

Li t r I Li 1 



n 



This isolates the UV divergence: when we use this small- s behavior in the integral (96) we 



get a divergence as e — > 0. To get a finite answer we just subtract off the contribution of Kq. 
This has the interpretation of renormalizing the various bulk and boundary cosmological 
constants. (Terms in Kq are proportional to the total volume, or the volumes of various 
walls or corners.) After making the subtraction we can set e = 0. So the renormalized 
answer is 



(3F 



1 



ds 



D 



where just for completeness 

K P (s,L) 

K D {a,L) 



L 



+ 



L 



— -V- 

47TS V 71 " 5 i 

71=1 



p vte D 



-L 2 n 2 /4s 



Li 



1 



'47TS 



(99) 



/47TS 2 Vtts ^ 

71=1 
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